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AFTER RUNAWAY: THE TRANS-HILL STAGE OF PLANETESIMAL GROWTH 

YORAM LlTHWICK^ 
Draft version March 19, 2013 

ABSTRAGT 

When planetesimals begin to grow by coagulation, they enter an epoch of runaway, during which the 
biggest bodies grow faster than ah the others. The questions of how runaway ends and what comes 
next have not been answered satisfactorily. Here we show that runaway is followed by a 'trans-hill 
stage' that commences once the bodies become trans-hill, i.e. once the Hill velocity of the bodies that 
dominate viscous stirring matches the random speed of the small bodies they accrete. Subsequently, 
the small bodies' speed grows in lockstep with the big bodies' sizes, such that the bodies remain in 
the trans-hill state. Trans-hill growth is crucial for determining the efficiency of growing big bodies, 
as well as their growth timescale and size spectrum. 

We work out the properties of trans-hill growth analytically and confirm these numerically. Trans- 
hill growth has two sub-stages. In the earlier one, coUisional cooling is irrelevant, in which case the 
efficiency of forming big bodies remains very low ( 0.1% in the Kuiper belt) and their mass spectrum 
is flat. This explains results from previous coagulation simulations for both the Kuiper belt and the 
asteroid belt. The second sub-stage commences when small bodies begin to collide with one another. 
CoUisional cooling then controls the evolution, in which case the efficiency of forming big bodies rises, 
and their size spectrum becomes more top-heavy. 

Trans-hill growth terminates in one of two ways, depending on parameters. First, mutual accretion 
of big bodies can become significant and conglomeration proceeds until half the total mass is converted 
into big bodies. This mode of growth may explain the size distributions of minor bodies in the Solar 
System, and is explored in forthcoming work. Second, if big bodies become separated by their Hill 
radius, oligarchy commences. This mode likely precedes the formation of fully-fledged planets. 



1. INTRODUCTION 

Our understanding of how planets form remains in- 
adequate. This hinders attempts to explain the many 
recent discoveries of extra-solar planets, proto-planetary 
disks and debris disks. Planet formation is often decom- 
posed into a number of stages. In the flrst stage, "plan- 
etesimals" form out of dust embedded in protoplanetary 
disks. How that happens is highly uncertain because it 
depends on complicated physics such as how particles 
stick and how they interact with a turbulent gas disk 
(see iChiang &: YoudinI l2010l for a review). The initial 
planetesimals are often assumed to be kilometer-sized, 
but they could be mu ch smaller or larger than that (e.g., 
iJohansen et al.l [20071 ). In the second stage, sometimes 
called coagulation, the planetesimals attract one another 
gravitationally, merge and grow. This stage is more eas- 
ily understood than the first, bec ause the dominant phys- 
ical process is simply gravity (see lGoldreich et al.ll2004R 
hereafter GLS, for a review). Nonetheless, past studies 
have given different, sometimes even conflicting, views 
on this process. What happens after coagulation de- 
pends on local conditions. Just beyond the snow-line, 
coagulation is thought to produce cores that then ac- 
crete massive gaseo us atmospheres to form gas giants 
([Pollack et al.l 119961 ) . In the terrestrial zone, coagula- 
tion produces dozens of sub-earth-sized bodies that then 
undergo a velocity instability once they have accreted 
half the planetesimals, lea ding to an epoch of large-scale 
chaos and giant impacts (jChambers fc Wetherilll 119981 : 
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iGoldreich et al.ir2004al ). And in the asteroid belt, Kuiper 
belt, and extra-solar debris disks, coagulation is thought 
to have been incomplete, perhaps because the coagulat- 
ing planetesimals were excited by exterior planets before 
forming planets themselves, or because the initial surface 
density was very low. 

It is the second stage, coagulation, that is the topic of 
this paper. Our goal is to build a theory that explains 
the properties of bodies that ultimately form, such as 
their number, size distribution, formation timescale, and 
efficiency, where the efficiency is the fraction of mass 
in the original planetesimals that ends up in big bodies. 
This theory can then be compared against observations, 
especially of the asteroid and Kuiper belts, whose large 
bodies are thought to be frozen remnants of the coagu- 
lation process. 

Coagulation is usually studied with numerical sim- 
ulations. Starting from an assumed initial state, the 
planetesimals begin to merge. The merging rate depends 
on gravitational focusing factors, which are functions 
of the relative speeds. The relative speeds are in turn 
affected by a variety of two-body processes such as 
viscous stirring, dynamical friction, and inelastic colli- 
sions. While there have been many s uch simulations, 
using a variety of techniques ( e.g.. |Greenberg et al.l 
1978[ IWetherifl fc Stewa^ 11981 IKokubo fc Idal 1 19981: 



WeidcnschiUing et al.' '199 7t iKenyon fc LuuT Il99£ 
Inaba et al. 2001: Morbidelli et al.l 120091 lOrmel et al.l 
2010bl : iWeidenschiUindWlH ISchlichting fc Saril l20Tll) 7 



the simulations are complicated, and hence it is often 
difficult to disentangle the various effects in order to 
understand the results and be confident that they 
are correct. Moreover, the results from different 
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groups do not always agree (e.g. IMorbidelli et~all[200l 
IWeidenschillind 120 lit) . Compounding the difficulty 
are a number of fundamental uncertainties — such as 
the unknown initial size and velocity distributions, 
how coUisional fragmentation occurs, and when are 
gaps opened in the circumstcllar disk. A theory for 
coagulation would be desirable to guide and interpret 
simulations, and to determine how sensitive the results 
are to assumptions. 

Early coagulation simulations led to the discov- 
ery of runaway growth, in which the size distribu- 
tion quickly develops a, tail extending to ext r emely 
large sizes (iSafronovl 119721 : iGreenberg et al.l 119781 : 
iWetherill fc Stewartl 119891 ^ ^Runaway occurs because, 
with gravitational focusing included, the growth rate 
of bodies {d\nR/dt) can be an increasing function of 
their radius R (see also Section [3.11 below). This im- 
plies that the largest bodies continue to double to infi- 
nite size before smaller ones double a single time. The 
properties of runa way g rowth have been studied ana- 
lytically Ce.g.. lLei[200a iMalvshkin fc Goodm"a3[200l . 
A variety of theories have been proposed to explain 
how runaway growt h ends, and what corr i es next (e.g., 
Ida & Makino Tl9M IKokubo fc Idal [T998t lOrmel et alJ 
2010k ,Schlichting fc Sarill2011[ ). It is often thought that 
runaway accretion is followed by self-regulated oligarchic 
growth, during which each big body heats its own food. 
We show here, however, that there is a critical interven- 
ing stage of growth, which we call the trans-hill stage. 

In the absence of an und ersta nding of the trans-hill 
stage, previous studies (^e.g.. lGLS[ ) could not explain two 
key results of conglomeration simulations: the size spec- 
trum and the formation efficiency of large bodies. Here 
we show that trans-hill growth is a critical stage for deter- 
mining these, and based on our analytical understanding 
we provide simple scalings to explain results of previous 
numerical simulations. 

During trans-hill growth, conglomeration can be cate- 
gorized as coUisionless or coUisional, with very different 
outcomes for formation efficiency and size spectrum. In 
the coUisionless case small bodies rarely collide, while 
in the coUisional one the small bodies' dispersion is 
reduced by frequent inelastic collisions. CoUisionless 
growth has been often simulated before, but without 
clear interpretation. CoUisional growt h, on the other 
hand , is just starting to b e explored (jWeidenschillingl 
1201 It [Shannon et al.l l2013b( ). Our study here provides 
a unified framework for both regimes. 

The structure of this paper is as follows. In Section [5] 
we present the equations of motion. We re-write these 
in Appendix A in a form suitable for numerical integra- 
tion or analytic solution. In Section [3] we review run- 
away growth, and show that it inexorably transitions into 
trans-hill growth. Wc then derive qualitatively the prop- 
erties of trans- hill growth. In Section 2] we present exact 
solutions of trans-hill growth, using numerical integra- 
tions as well as analytical self-similar solutions (derived 
in Appendix B). Section [5] describes what comes after 
trans-hill growth. Section |6] examines our assumptions 
and delineates their range of validity. Section [7] discusses 
applications of our theory to the asteroid and Kuiper 
belts, and compares with a number of earlier papers. 
Section [8] svmimarizcs. 

.Shannon et al., (,2013allb[ ) present particlc-in-box simu- 



lations in both the coUisionless and coUisional regimes, 
without the restrictive assumptions made in the present 
paper. The results there confirm and refine those of this 
study. The coUisional paper in particular focuses on the 
formation of the Cold Classical Kuiper belt. It presents 
a new picture where KBOs form out of a very low mass 
planetesimal disk (the "minimum mass Kuiper belt"), 
100 times less massive than the minimum mass solar neb- 
ula. 

2. ASSUMPTIONS AND EQUATIONS OF MOTION 

Wc examine first the interactions between two groups 
of bodies, big ones and small ones, before proceeding to 
consider a distribution of bodies. T his 't wo-groups' ap- 
proximation has been described by IGLSI and we follow 
their notation. Big bodies have radius R and surface 
density S; small bodies have radius s, surface density ct, 
and velocity dispersion u. All bodies have the same bulk 
density, with p ^ 1 g cm~'^. We make the following as- 
sumptions, and then check for self-consistency in Section 

El 

1. The small bodies' random speed u satisfies 

< u/vH < a-'^^ , (1) 

where a is the ratio of all bodies' physical radius to 
their Hill radius {R/Rh) and vh is the big bodies' 
Hill velocity. Explicitly, 



where the expression Rq / a (the ratio of the Sun's 
radius to the semimajor axis) follows from the fact 
that the Sun's density is comparable to that of solid 
bodies. In the asteroid belt a ^ 2 x 10~^, while in 
the Kuiper belt a ^ 10^"*. In addition, 

Vh ^ R\/ Gpa , (3) 

The small bodies' speed can be sub- or super- 
hill with respect to the big bodies [u < vh or 
u > Vh), subject to the constraints that u is less 
than the surface escape speed from the big bod- 
ies {u < a^^/^vn) and that small bodies are not 
accreted in the "very thin disk" accretion regime 
(m > a^/^w//). 

2. The big bodies' random speed v is sub-hill {v < 
Vh)- 

3. S <C cr, and so a is essentially the total surface den- 
sity, and is treated as a constant during the growth. 
More stringently, we assume that big bodies grow 
only by accreting small ones. This is certainly true 
at early times, when there are only a few big bod- 
ies. 

4. The small bodies' size s is a constant parameter. 
Small bodies do not grow, and they also do not 
fragment in collisions. 

5. More than a single big body dominates viscous stir- 
ring at a given distance from the star, and hence 
the big body number distribution can be treated as 
a continuous function. This assumption is violated 
when oligarchy commences. 
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With the above assumption s, big bodies grow by ac- 



creting smah ones at the rate ()GLi 



1 dR 



an 



vh/u, 



if u > Vf{ 
if u < vh 



(4) 



where f2 is the orbital angular speed around the Sun. 

For the evolution of u, small bodies are damped by 
inelastic collisions amongst themselves, an d are v iscously 
stirred by the big bodies. The net rate is (|GLSD . 



u dt 



ps pR 



vh/u, 



if u > Vh 
if u < Vh 



(5) 



If the first term on the right hand side is important the 
system is coUisional, otherwise it is coUisionless. 

The above equations for two groups of bodies are easily 
extended to a continuous distribution of big bodies. We 
denote the cumulative number distribution of big bodies 
as N{> R). Because we neglect accretion of big bodies, 
the number of big bodies remains constant as they grow. 
In other words, N satisfies the continuity equation 



dt ^ dR ' 



(6) 



where V = dR/dt is given in Equation (j4|). In place of 
Equation ([5]), we replace S dS, then set 



d{J:/a) 
d\nR 



= R 



d{N/v) 
dhiR 



(7) 



and integrate over dR. Here 77 is a constant; we shall not 
need its value because it is only the ratio N/rj that is 
dynamically significant. 

Equations (jl])-!?]) are the equations of motion. We 
solve them both analytically (Appendix B) and numeri- 
cally (see method in Appendix A). But before presenting 
the exact solutions, we first use simple analytical argu- 
ments to show that runaway growth inexorably converges 
towards 'trans-hill growth'. This realization gives rise to 
a number of our main results on the mass spectrum and 
formation efficiency (Section [3]), as verified by our exact 
solutions (Section |4]). 

3. FROM RUNAWAY TO TRANS-HILL GROWTH 
3.1. Runaway 

We summarize the tr aditional picture of h ow runaway 
growth proceeds (e.g., Ilda &: Makinol [l993l ) . It is typ- 
ically assumed that bodies of some characteristic size 
(here, s) emerge from a dissipating protoplanetary gas 
disk. The value of s is highly uncertain, since it is not 
even understood how the bodies formed. 

On the timescale that bodies of size s collide, ps/ (an), 
they stir each other up to their surface escape speed, 
u ^ s^/Gp ~ a~^/'^VH\s (Equation (O), where vh\s 
is the Hill velocity for bodies of size s. On the same 
timescale, they grow in mass by accreting each otherQ 
The evolution of those bigger bodies (radius labelled by 
R) proceeds in the runaway regime. This is because the 
growth rate when u > vh is an increasing function of 

^ It is possible that some bodies grow to very large sizes before 
u is stirred up to the small bodies' escape speed, depending on the 
details of how the small bodies formed. We discuss the effect of 
this on trans-hill growth below. 



R, i.e. dlnR/dt cx R (Equations © and Thus, 
comparing two big bodies with different values of R, the 
bigger one will double in size faster than the smaller 
one, and then will double again even faster. It will ulti- 
mately reach infinite size — or violate one of the assump- 
tions made — before the smaller one has doubled a single 
time. 

3.2. From Runaway to Trans-hill Growth 

Runaway growth requires u > vh- Since vh increases 
linearly with i?, when runaway bodies get big enough 
their Hill velocity can be sufficiently large that u < vh- 
In that sub- h ill case, big bodies grow in the "neutral" 
regime i.e., since the growth rate is independent 

of R (Equation Q), the distribution function of big bod- 
ies N{> R) maintains its shape while moving to larger 
R. 

In general, some big bodies grow in the runaway 
regime, and others in the neutral regime, depending on 
their radius: 



runaway growth, if i? < i?tian 
neutral growth, if i? > -Rtran , 



inhere the trans-hill radius is 

Rtmn — 



y/Upa 



(8) 
(9) 

(10) 



which is the radius of the big bodies that have Hill ve- 
locity equal to the small bodies' speed {vH\Rtr^„ = "")• 

As we now show, the big bodies' distribution func- 
tion is always driven to the trans-hill state, defined as 
i?stir ~ ^tran, where i?stir IS the radius of the big bod- 
ies that dominate viscous stirring. Subsequently, the 
small body velocity dispersion and the size of the bod- 
ies that dominate stirring grow in lockstep, maintaining 
i?stir Rtrun- Our argument proceeds by considering the 
two initial situations, i?stir < ^tran or i?stir > i?tran, in 
turn. 

• If initially i?stir ^ ^tran (Fig. [U, as is relevant 
for the early phases of planetesimal growth, the 
small bodies are super-hill with respect to the stir- 
rers. All big bodies with R > i?stir repeatedly dou- 
ble in size before the stirring bodies (i?stir) have 
doubled onceO This runaway produces new bodies 
that contribute more to stirring (cx Si?'^, Equation 
([5])) than the original stirrers and these new bodies 
come to dominate the stirring. Although i?tran also 
grows due to viscous stirring, i?stir grows faster, as 
we presently show. The growth of i?stir is at least as 
fast as the size doubling rate at i?stir , d In i?stir /dt > 
d\nR/dt\i^-ii^^^^. Meanwhile, the heating of small 
bodies increases i?tian as din Rtj-i^^/dt ^ dlnu/dt 
(Equation ([5])). We first focus on the coUisionless 
case, inserting s = 00 into Equation ([5]) to get 



din Rstir/dt 

din Rtran/ dt 



> 



Rtia 
Rsti 



where Sgtir is the mass density at i?stir- If the ef- 
ficiency of formation is limited by Sstir/f ~ ct, as 

^ Even though very large bodies, with radius above -Rtran, grow 
in the neutral regime, their doubling rate also exceeds that of the 
stirring bodies. 
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Figure 1. From Runaway to Trans- hill Growth. Result of a nu- 
merical integration of the coUisionless equations of motion (Equa- 
tions Q^lO with s = oo and a = 10^^). The leftmost panels 
(top and bottom) show the initial conditions, where the bottom 
panel is the big body distribution function, and the top panel is 
the stirring rate per In R (the right-hand side of Equation l(5} with 
S — >■ \dT,/dln R\). The red cross marks Rstin defined to be the 
maximum of the curve in the top panel. The blue circle marks 
R tran = V GpOL. Initially u is the escape speed from bodies 
with R = 1 (i.e. -Rtran = a"^''^ = 100). The panels at subsequent 
times (r = 70, 500, 10*^, with t the scaled time, t = taQ/{pa); see 
Equation l lAll l) show that the system evolves from -Rstir ^Jtran 
to the trans-hill state Ratii ~ ^?tran, and thereafter remains trans- 
hill. 

we show below to be true for the colhsionless case, 
with time i?stir catches up to i?tran- An alternative 
argument showing that i?stir grows faster is that if 
there were no bodies bigger than the original stir- 
ring bodies, then, as long as the small bodies are 
colhsionless, the stirring bodies grow at the same 
rate as u, i.e. dlnRstii/dt ~ dlniitran/c^^- But by 
the nature of runaway growth, the bodies bigger 
than the original i?stir must grow faster than that. 
We conclude that, in the colhsionless regime, the 
ratio i?stir /-Rtran approachcs unity. The simulation 
in Figure [1] confirms this. 

The same conclusion applies to the collisional case 
as well, because collisional damping forces u (and 
hence i?tran) to grow even more slowly than in the 
colhsionless case. Hence i?stir catches up with -Rtran 
even more quickly. 

• If initially i?tran ^ -Rstir, the stirring bodies grow 
in the neutral regime. Focusing first on the colli- 
sional case. Equation ([Sj with u < vh implies that 
Sstir oc M (X -Rtran- Bccausc the Stirring bodies grow 
neutrally, the shape of the distribution function 
near -Rgtir does not change as -Rstir grows, and hence 
the number of stirring bodies does not change, i.e., 

-Rstir- As a resuh, -Rstir/-Rtran OC 1/-Rstir: 

which decreases with time towards unity. The same 



conclusion applies to the colhsionless regime, be- 
cause in that case u, and hence -Rtran, will increase 
even faster. 

3.3. How Trans-hill Growth Proceeds 

We have shown that during runaway the big bodies' 
distribution is driven towards the trans-hill state 



Rs 



Rt 



or equivalently 



VH 



(11) 



(12) 



and that after that the big bodies continue to grow in 
this state. We call this new stage of growth trans-hill 
growth. 

During trans-hill growth, all bodies with R < -Rstir 
do not grow significantly because those with R > -Rstir 
double faster than they do, i.e. the distribution becomes 
frozen in time at small R. Furthermore, bodies with 
R > -Rstir follow neutral growth, and their distribution 
maintains it shape. Throughout, -Rtran ^ -Rstir, and both 
increase with time, following (Equation |4]) 



aft 

-Rstir ^ 1 

pa 



(13) 



We may also evaluate the mass fraction in bodies that 
dominate stirring, Sgtir/c, as a function of -Rgtir- We 
find, dropping the "stir" subscripts. 



a 



a, coUisionless (-R < s/a) 

a^R/s, collisional (-R > s/a) 



(14) 



where the colhsionless expression above follows from set- 
ting dlnu/dt ^ d\n Rstir /dt and dropping the collisional 
term in Equation ([S]); and the collisional expression fol- 
lows from dropping the left-hand side of Equation (IS|)Q 
The transition from colhsionless to collisional trans-hill 
evolution occurs once -Rstir exceeds a critical size 



i?stlr>i?iTr^=s/a 



(15) 



because beyond this size, small bodies coUisionally cool 
faster than the growth of trans-hill bodies. This transi- 
tion occurs at t ~ ps/{aVl), i.e., the small bodies' colli- 
sion time. 

Following conventional practice, we define the power- 
law index q of the differential number distribution 



dR 



cxl](> R)R- 



in which case we have 



4, colhsionless (-R < s/a) 

3, collisional {R> s/a) 

4. EXACT SOLUTIONS 



(16) 



(17) 



We present exact solutions of the equations of motion. 
We first summarize the analytical self-similar solutions to 

* Because the spectrum is frozen at small R, Equation I I14I I gives 
not only the temporal evolution of Sstir, but also the frozen spec- 
trum at -R < -Rstir — that is why we drop the "stir" subscript. 
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these equations (derived in Appendix B), and then inte- 
grate the equations numericahy for a number of cases — 
both when the initial conditions are as given by the self- 
similar solution, and for more realistic initial conditions. 

4.1. Analytical self-similar solutions 

It is shown in the appendices that the equations of 
moti on (Equations (H])-©, or equivalently Equations 
(|A1|) - (|A4|) ) admit self-similar solutions when the radii R 
are scaled relative to i?tran, and when the stirring is either 
in the collisionless or collisional regime. In either case, 
there is a single free parameter, the power-law exponent 
for bodies that grow neutrally, i.e., the value of 7 such 
that N oc R~'' at i? > -Rtran- The value of 7 remains 
frozen in time by the nature of neutral growth, as long as 
the assumptions of Section [2] remain valid. However, the 
value of 7 in real disks is difficult to ascertain from first 
principles, and likely depends both on how the bodies 
form and on the early stages of coagulation. Fortunately, 
the solutions are quite insensitive to 7 as long as 7 ^ 3, 
or equivalently as long as stirring is dominated by bodies 
at finite sizes, not by those with R = 00. 

The exact sel f-simi lar solutions (displayed in Equa- 
tions p8)) and (|B12p and graphed in Figures [2] and O 
confirm what was derived qualitatively in Section [3] In 
particular, from the functional dependence of N on R, 
one sees that stirring is maximized near R i?tran, and 
hence i?stir ^ ^tran- In addition, the differential mass 
fraction at R^ Rtmn is consistent with Equation (fT4| : 



dluRa 
d S 



-3gnocoi(7)a, colhsionlcss (18) 



R 



dlnR (7 



= —2gcoi{'y)a — , colhsional 



(19) 



The order- unity functions (7nocoi(7) and .gcoi(7) are de- 
fined in Appendix B. They asymptote to (?nocoi(7 — >■ 
00) = .gcoi(7 ^ 00) = 10. 

4.2. Numerical solutions 

Figure [2] shows a numerical integration of the equations 
of motion in the collisionless case (s = 00), initialized 
with the collisionless self-similar solution (Equation IB8|) 
at scaled time r = 1. It is apparent from the figure that 
the evolution remains self-similar at all times, and agrees 
with the analytic expression (Equation (jBSI) ). The mass 
distribution per logarithmic bin is indeed constant, in 
agreement with Equation p7p . 

Figure |3] shows an integration for th e collisional case. 
It is initialized with Equation (jB12[) at r = 1, has 
s = 0.1a X -Rtran|r=i, and other parameters as before. 
The numerical integration agrees with the analytic ex- 
pression, and shows that the mass distribution is top- 
heavy with dN/dR oc R^^ (Equation p7| ). 

In Figure |4l we experiment with a more commonly 
adopted initial spectrum — one that is more strongly 
peaked at a single size (R = 1). We perform two in- 



tegrations, both of which have a = 10 



1 (i.e., s 



is the same as the location of the peak in the initial size 
distribution), and initially N oc at i? > 1. The two 
integrations differ in the initial u: one has initial veloc- 
ity equal to the escape speed of bodies with R = 1, and 
the other has u ten times lower. Over time, both con- 
verge to the collisional self-similar solution. The low u 
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Figure 2. Collisionless evolution of Equations (|4]|-([71l, initial- 
ized with the self-similar size distribution at scaled time t = 
taQ/{pa) = 1; -R is measured in units of the initial value of -Rtran- 
Top panel: The cumulative number distribution of big bodies (nor- 
malized by constant rj) is plotted at three subsequent times. The 
numerical and analytic self-similar solutions agree. The initial pa- 
rameters are s = 00, a = 10^'*, and 7 = 7, i.e., oc R^^ at large 
_R. Bottom panel: The differential mass distribution, derived from 
the top panel via Equation {TJ. At each time, the radius of 
bodies that dominate stirring (Rstir) is near the break in A'' (or S) 
from one power-law slope to the other, and the small body speed 
is such that u is trans-hill relative to those stirring bodies (i.e., 
^?tran/i?atir = order-unity constant). Collisionless evolution leads 
to a size spectrum dN/dR oc i?-* at sizes below -Rgtin and the 
efficiency of forming large bodies is low, d{^/a)/dlnR ~ 10a. 

case undergoes a temporary phase of collisionless trans- 
hill growth, when the big bodies have yet to grow to 
R > s/a. By contrast, the high u case takes longer 
to reach the trans-hill solution because its initial condi- 
tion is more discrepant from trans-hill. In fact, this case 
completely skips the collisionless trans-hill regime. We 
have also experimented with different values for the ini- 
tial power-law slope of N. The main resulting difference 
is the value of R at which the collisionless trans-hill solu- 
tion commences. For example, when N oc R^^ aX R > 1 
(rather that R~^ as in Figure[3]), then the high u solution 
follows the collisionless theory from R = lO^'^-lO"*. 

In summary, as long as the evolution is collisionless, 
the efficiency of forming big bodies is small, with S/ct ~ 
lOa, independent of R. But once the evolution becomes 
collisional (i?stir > s/a), the efficiency grows towards 
unity as i?stir increases. 

5. AFTER TRANS-HILL GROWTH 

5.1. Growth by accreting big bodies: Equal accretion 

One of the ways trans-hill growth can end is when ac- 
cretion of big bodies becomes important, violating as- 
sumption 3 in Section [21 To establish when that occurs, 
wc calculate the accretion rate of big bodies, which re- 
quires knowing the big bodies' random speed v. Balanc- 
ing dynamical friction damping due to small bodies with 
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Figure 3. CoUisional evolution: Similar to Figure|2] but showing 
the coUisional case with parameters s = 0.1a (in units where the 
initial -Rtran = 1) and, as before, a = 10~^ and 7 = 7. The size 
spectrum is dN/dR <x R^'^ at sizes below -Rstin a^nd the efficiency 
of forming large bodies is much higher than the coUisionless case. 



viscous stirring due to other big bodies yields (jGLSD 

1 dv n 



dt pR 







(20) 



Therefore v/vh ~ S/cr^ 1. T he growth rate by accret- 
ing big bodies with v < vh is (jRafikovll2003al GLS) 



1 dR. 
R~dt^^'e 



'pr' 



'3/2 



(21) 



which is larger than the usual (isotropic) sub-hill accre- 
tion formula (e.g., Equation (j4|)) because big bodies lie in 
a zero-inclination disk. Comparing growth by accreting 
big bodies with that by accreting small bodies, 



c? In i?/c?t I big S 

In I small cr 



a 



-1/2 



1/2 



a^/^Rstn-/s, 



coUisionless 
coUisional 



(22) 
(23) 



using the trans-hill expressions. We conclude that in the 
coUisionless regime it is always safe to ignore accretion 
by big bodies, in disagreement with iSchlichting fc Sari 
()201ll) (see also Section [71). However, in the coUisional 
regime, big body accretion becomes important once i?stir 
exceeds a critical value 



^ T-,(b.b. accrete) 



stir 



v-3/2 



2800 km ( 



V 1 km 



VlAU 



3/2 



(24) 



This occurs at time t ^ a ^/"^psj (crf^), and at that time 
the fraction of mass in big bodies (i.e., the efficiency) is 
E/a ^ ~ 7%(a/lAU)-i/2. 



10° 



10-1 



K 10- 



^ 10-3 



10- 



iiij 1 1 null 1 1 iiiiiij 1 1 iiiiiij 1 


Illlllj 1 1 llllll 


1 1 Illlllj 1 1 Illlllj 1 1 Illlllj 1 Mill 






^ / — 

to ■ 

. o, • 


\\\ low u{t 




UlllolUlllCba 1 


— \\\ 1 




theory \ 


\ \high u( 
_ 'r = 0\\ 


/=0) 


s/a 1 


ml 1 1 1 1 1 


1 nI 


^ 1 1 1 Mini 



10-^10° 101 104 105 106 107 

R 

Figure 4. Prom runaway to trans-hill. Two integrations were 
initialized with nearly all bodies at the same size, R= \ (the black 
profile marked r = 0). Numerical integrations of the equations of 
motion are plotted at scaled time t = 2.5 X 10®. For the red (high 
u) curve, the initial u is the surface escape speed from bodies with 
R = \ (i.e., -Rtran = ^ / = 100); for the blue curve, it is ten 
times smaller (-Rtran = 10). Both integrations begin in the run- 
away stage, an d converge onto the coUisional self-simil ar s olution 
(Equation JTOl l. dotted line) when R> s/a (Equatio n ([ril ). The 
low u case follows the coUisionless solution (Equation 11181 . dashed 
line) for some time before coUisional cooling becomes important. 
This explains the fiat mass spectrum from .R = lO'^ — 10*. The 
high u case takes longer to reach the self-similar solution because 
its initial conditions differ more from the trans-hill state; in fact, 
it effectively skips the coUisionless solution entirely. 



As shown in [Shannon et al.l ()2Q13bD , after coUisional 
trans-hill accretion ends, an epoch of "equal accretion" 
begins, during which big bodies grow by accreting com- 
parable mass in big and in small bodies. Equal accretion 
terminates when half of the mass has been converted into 
big bodies, E ~ cr. 

The scenario outlined above is applicable as long as 
oligarchy has not yet begun. We discuss oligarchy next. 

5.2. Oligarchic Growth 

As accretion proceeds, the number of stirring bodies 
decreases, and hence neighboring stirrers become increas- 
ingly separated. Eventually, they become so separated 
that each small body is predominantly stirred by, and 
accreted onto, a single big body. When that happens, 
assumpt ion 5 is violated (Sec tion [2l) and oligarchy com- 
mences (IKokubo fc Idal [l998l . GLS). In oligarchy, the 
nature of growth is modified. Instead of a continuous 
size spectrum, the largest body in each radial annulus 
separates from the size spectrum of smaller bodies. The 
velocity dispersion u is no longer trans-hill. 

The value of i?stir at which oligarchy starts depends 
on the size spectrum, for which we now have a simple 
model. For trans-hill velocity dispersion, oligarchy be- 
gins when the separation between adjacent big bodies 
(Aa ~ pR'^/ (aS)) exceeds their Hill radius {Rh ~ R/a). 
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0P = 



If we define the oHgarchy parameter 

Aa 

_ paR^ 

then oUgarchy begins when 0P~ 1. Inserting the efii- 
ciency of big body formation (Equation ([T4| ) into the 
above, we find that trans-hill growth transitions to oli- 
garchy once i?stir exceeds 



(25) 



1. The requirement that u satisfies Equation ([1} is 
equivalent to insisting that we only consider big 
bodies with a^/^ < i?/i?tran < Bodies that 

violate these restricti ons grow more slowly than 
Equation (g]) predicts (jGoldreich et al.ll2004'b| ). but 
this has little effect on trans-hill growth. 



The assumption that v/vu 
Equation 



R. 



(oligarchy) 



150km. 



f) ^« 120km. 



.1/2 



collisionless 
(26) 

colhsional 4. 



< 1 is confirmed by 



We ignored growth by accretion of big bodies. 
When that assumption is violated, equal accretion 
begins fSection lS.ip . 



where (Tig = a/{lQg/ cm^). For the MMSN density pro- 
file, CTie ~ (a/lAU)^'^/^ with an enhancement beyond 
the snow line by a factor of ~ 5. The collisionless ex- 



pression above applies as long R^^^^^'^'^^^'^ 



<R 



(col) 



s/a; 



otherwise the collisional one applies. However, we cau- 
tion that depending on parameter values the equal ac- 
cretion stage discussed in Section [O] may begin earlier 
than oligarchy, which would lead to a different transition 
criterion. 

Once oligarchy takes hold, neighboring regions evolve 
independently under the stirring of t heir r espective oli- 
garchs, as discussed in Sections 9-10 of lGLSI . Neighboring 
oligarchs likely converge in size, their battling leading to 
scattering or merging. The evolution depends to a large 
degree on the fate of the small bodies. At late times, 
small bodies collide at such high speeds that they almost 
certainly fragment. As the small bodies grind down each 
other, the evolution can become increasingly collisional. 
The small bodies may cool and the oligarchs may carve 
gaps around themselves, effectively sabotagi ng their ac- 
cretion pla fikov 2003bi: iLevison et al.|[2010l ). We defer 
considerations of these dynamics to future work. 

If the oligarchs eventually reach the isolation mass, 
they may undergo orbital instability and experience gi- 
ant impacts. The timescale to form planets of size R 
is 



t 



impact ■ 



pR 



10» 



R 

R(n 



lAU 



3/2 



yrs, (27) 



which is very long in the outer solar system. 

Our criterion for o ligarchy (0P~ 1) differs from that of 
lOrmel et al.l ()2010aD . They stipulate, based on empirical 
evidence from Monte Carlo simulations, that oligarchy 
begins when the stirring rate by one single large body 
equals its growth rate, and that, in turn equals the small 
body collision rate. Their definition of oligarchy differs 
from ours: their oligarchy begins when the ratio of small 
body random velocity to Hill velocity of the biggest body 
reaches a minimum. In fact, their oligarchy is similar to 
(but still somewhat different from) the trans-hill phase. 
To some extent, this is a matter of how one defines oli- 
garchy. But irrespective of definitions, there should be 
a transition from trans-hill to oligarchic behavior when 
Aa - Rh- 

6. EXAMINING ASSUMPTIONS 

We examine the assumptions made in deriving the 
properties of trans- hill growth (Section [2]): 



We assumed that s is constant. In truth, at late 
times u can become sufficiently large that collisions 
between small bodies fragment them. We do not 
treat this in detail because the physics of fragmen- 
tation is complex. We note, however, that since 
bodies of smaller sizes are typically more resistant 
to fragmentation (for s ^ 1 km), the result of frag- 
mentation is likely that s is a decreasing function of 
the stirring bodies' sizes. This will alter the spec- 
trum of trans- hill accretion (Equation (|14p ) in that 
one should replace the s in that equation with the 
function s(i?stii). A similar remark applies to other 
mechanisms for damping u, such as gas drag. Since 
s only appears in the equations of motion through 
its damping effect on u, gas drag may also be mod- 
elled by adopting a much smaller effective s. 

5. We ignored oligarchy. The transition to oligarchic 
stage is described in Section [5^2] 

7. APPLICATIONS AND COMPARISONS 
7.1. The Asteroid Belt 

Here we compare our results with published coagula- 
tion simulations for the asteroid belt. At a distance of 2 
AU, a ss 2.5 X 10"'^, so the transition from collisionless 
to collisional evolution occurs when large bodies grow 
beyond s/a ^ 400 km(s/l km). 

The uppermost end of the asteroid mass distribution 
is roughly flat, i.e. dT,/d\nR ~ const for 100km < 
R < 300km But for R 5, 100km, £ falls off with de- 
creasing R (jJedicke et al.ll2002t iBottke ej al. 200 51, tak- 
ing a constant alb edo of . 04). Both iMorbidelli et al.l 
(|2009( ) and Weiden schillinel (|2011| ) attempt to reproduce 
the observed distribution, including the 100km "bump," 
by running particle- in-a-box coagulation simulations that 
are initialized with single-sized small pl anetesimals. But 
the tw o reach opposite conclusions. IMorbidelli et al.l 
(|2009f l. with small plantesimal sizes ranging from s = 0.6 
to 6 km, fail to produce the 100km bump: the final dis- 
tributions they find are roughly flat (q ^ 4) for a range 
in R that extends an order of magnitude below 100km. 
But Weidenschilling succeeds when using s « 0.1 km. He 
attributes the difference to the erroneous neglect of sub- 
hill accretion by Morbidelli et al. In the following, we 
focus on explaining Weidenschilling's results. 

^ An additional diff erence between the simulations of 
IMorbidelli et al.l | |2009|'| and lWeide nschillini 120111) is that the for- 
mer initialize velocities to be the Hill velocity of the initial bodies, 
whereas the latter initialize them to be the escape speed. This 
might tend to flatten out the post-runaway bump in the Simula- 
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Weidenschilling's bump at lOOkm appears to be due 
to the transition from runaway to collisional trans-hill 
growth. These simulations resemble the red curve in 
Fig. |4] — with a transition directly from runaway to col- 
lisional trans-hill — because the initial velocity is the es- 
cape speed from seeds. And like the red curve that shows 
a transition around R ^ s/a, the final size distribution in 
iWeidenschillingj (|2011[ ) also shows a bump around 50km 
(his Fig. 8), produced when the small body dispersion 
goes from super-hill (runaway) to trans-hill. The value 
of S/cr at the location of his bump, and the slope to the 
right of the bump, are all roughly consistent with what 
we derive in this paper. 

If the observed 100km bump in the asteroid belt is in- 
deed produced by the transition from runaway to trans- 
hill growth, that would imply that the initial mass of 
small plantesimals, relative to that in current big bodies, 
was at least tr/S ~ l/(10a) ~ 50 times greater (Equation 
(jl8p ) — or even greater if some a steroids were dynamically 
ejected (jMorbidelli et al.ll2009() . It would also imply that 
the initial small bodies had radii ~ a x 100km- O.2kn0 
One interesting implication of such a scenario is that the 
vast majority of these small bodies could not be dynami- 
cally ejected, because if they were that would correspond- 
ingly reduce the number of big bodies too. Instead, they 
had to be ground down to dust which was then elimi- 
nated by radiation forces. Whether this can occur has 
yet to be examined carefully. But we note that by the 
time the asteroids grew to — 100km, the small bodies 
would have begun to collide with one another with col- 
lision speeds ?^ff|_R=iookm ~ 5m/sec. Such speeds might 
have been sufficient to grind all the planetesimals and 
eliminate them, preventing further growth of asteroids. 

7.2. The Kuiper Belt 

With a value of a « 10^^ for the Kuiper belt, 
the transition to trans-hill collisional growth occurs at 
10"* km(s/l km). So with the typical choice s = 1km 
and the largest bodies observed at 1000 km, there is lit- 
tle wonder that simulatio ns to date have not pro bed the 
c ollisional regime (but s ee [Shannon et al.ll2013bf ). 

iKenvon fc Luul ()1999D perform particle-in-box simula- 
tions for the formation of Kuiper belt objects. Starting 
from bodies of size s = 80 m, their simulations pro- 
duce big bodies with a flat mass distribution {q = 4) 
for 3km < R < 1000km (see, e.g. their Fig. 8). This 
is roughly consistent with the observed m ass distribu- 
tion in the Kuiper belt at large sizes (e.g.. iFraser et all 
[2OT0t ). and hence may be taken as evidence that the for- 
mation of Kuiper belt objects has been solved. However, 
recent data have cast doubt on such a simple picture 
(jShankman et al.ll201^ . 

Our theory explains some aspects of Kenyon & Luu's 
simulations. The growth they witness should entirely 
be in the coUisionless regime, and the mass spectrum 
should be flat, as indeed they observe. However, the 
efficiency that they find is — 1%, whereas we predict 
— lOa — 0.1%. We have been unable to resolve this 

tions of MorbidcUi ct al., similar to how, in Figurc|4l the blue curve 
is flatter than the red. 

® The initial small body radii would be less than 0.2 km by a 
factor of ~ 5 if the MMSN budg et of gas was still present when 
the bump was formed (e.g.. ITJEH) . 



discrepancy. In [Shannon et all (l2013al ). we repeat their 
si mulations and find an effi ciency of '--^ 0.1%. 

iSchlichting fc Saril ()201lD consider coUisionless coag- 
ulation both analytically and with numerical simula- 
tions. They conclude that the flat mass spectrum ob- 
tained in colUisionless simulations arises because big 
bodies grow equally by accreting big and small bodies 
("equal accretion"). By contrast, we have shown that 
it is due to trans-hill growth. In fact. Equation (p3|) 
shows that trans-hill growth in the coUisionless regime 
prevents equal accretion from occuring. Nonetheless, the 
efficiency they derive {a^^^), happens to be close in nu- 
merical value t o our efficiency of 10 a for Kuiper belt 
parameters. In [Shannon et alI(|2013aD . we confirm with 
numerical simulations that account for order-unity co- 
efficients that it is trans-hill growth rather than equal 
accretion that sets the size spectrum in the coUisionless 
regime. 

A troubling concern with low-efficiency coUisionless 
growth of KBOs is that one needs to grind down 99.9% 
of the small bodies to dust and blow it out with radia- 
tion pressure. As we show in an upcoming paper, this 
is almost impossible to do. This concern is alleviated if 
the initial planetesimal size is significantly smaller than 
1 km. The formation efficienc y can then be significa ntly 
boosted by collisional growth ([Shannon et al.|[2013b[ ). 

7.3. \Ormel et all miOW ) 

lOrmel et al.l (|2010br ) perform a comprehensive study 
of conglomeration, at semimajor axes ranging from 1 
AU to 35AU. Their numerical algorithm allows them 
to follow the evolution from runaway through oligarchy. 
We observe that there is a distinct trans-hill phase in 
their simulations (e.g.. Figs. 8 & 11 of that paper). 
With their parameters of s = 7.5km, a — 16g/cm^, 
one expects the transhill evolution (coUisionless below 
-Rstir = 1500(a/lAU) km) to lead to a characteristic 
size spectrum of 9 = 4 and a formation efficiency of 
S/cr — 0.05(a/lAU). These are indeed observed in their 
resultsQ Above R = 1500 km(a/lAU), the evolution is 
collisional. For some of their simulations (the ones at 
small a), oligarchy enters at around the same point. So 
instead of observing a continuous mass spectrum of g = 3 
(Equation ([T7)) ). they find that a single body that grows 
to large sizes. 

8. SUMMARY 

Runaway growth ends once the stirring bodies become 
trans- hill (f/zlstir ~ m, or equivalently i?stir ^ uj ^JGpct). 
Afterwards, i?stir and u grow in unison, and the stir- 
ring bodies remain trans-hill. For i?stir < s/a, trans-hill 
growth is coUisionless. The mass spectrum is fiat (g = 4), 
and the efficiency low, '-^ 10a. But for i?stir > s/a, the 
evolution is collisional and the efficiency grows in propor- 
tion to i?stir, with a size spectrum that has q = 'i. The 
time to reach the collisional transition is comparable to 
the small bodies' collision time {ps/a^l). Our numeri- 
cal simulations and self-similar calculations confirm these 
results, subject to the assumptions made in Section [2j 

^ Although lOrmel et al.l Il20l0bl1 fit their result with a power-law 
slope of q = 5.5 (p = — 2.5 in their notation), it appears that the 
big bodies are better characterized by a flat mass distribution of 
q = 4 (see, e.g., their Fig. 13). 
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The simulations also show that where the collisionlcss 
trans-hill spectrum begins depends on the unknown ini- 
tial conditions — especially the initial size spectrum and 
It. But the coUisional trans-hill spectrum invariably be- 
gins once i?stir > s/a. 

Having described trans-hill growth, we briefly dis- 
cussed what comes next, when one or more of the as- 
sumptions in Section [2] breaks. CoUisional trans-hill 
growth ends once i?stir ^ sa~^^^ (efficiency E/cr > \/a), 
because once that happens mutual accretion of big bod- 
ies becomes import ant (Section [5t. As w e show in a 
forthcoming paper jSh annon eFaLl [20130 ). this regime 
is characterized by equal accretion of big and small bod- 
ies, and can lead to order-unity efhciency. However, the 
mass spectrum below i?stir is no longer frozen, and hence 
if this regime is reached it can wipe out the mass spec- 
trum laid down during trans-hill growth. 

Trans- hill growth can also be terminated when big bod- 
ies become separated by a Hill radius, in which case oli- 
garchy sets in. Another complication is coUisional frag- 
mentation. We leave more detailed consideration of these 
two effects to future work. On the other hand, the com- 
plication of gas damping is easily incorporated by using 
an effective s. 

Our theory explains results from previous studies of 
the asteroid and Kuiper belt. Most simulations to date 
adopt parameters relevant for collisionlcss growth, typ- 
ically with initial seeds s ^ 1km. For the Kuiper belt, 
this means the growth is collisionlcss at all times and 
the efficiency is limited to lOa ~ 0.1%. So if the Kuiper 
belt has been severely dynamically depleted, it would 
initially have had to contain multiple times the MMSN 
mass. Alternatively, the efficiency of formation can ap- 
proach unity if the initial seed size is small. Wc explore 
how this may explain the formation of Kuiper belt in an 
upcoming publication. 
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APPENDIX 

A. NUMERICAL SOLUTION OF THE EQUATIONS OF MOTION 

In order to numerically integrate the equations of motion (Equations (|4])-(l7])), we first re- write the equations in 
terms of the trans-hill radius i?tran = uj^fGpa. (Equation (jlOp) and the rescaled time, 



T = 1 

pa. 



which has units of length. Equation (|4|) becomes 



Cfi? ^ f (i?/i?tran)', if i? < i?tra„ 
dr 1 RIR if i? > i?tran 

and Equation ([Sj , with the replacement of Equation ([7]) , becomes 

1 di?tran _ a 1 / f dN f dN 

iJtran dr S a Vfl<flt„„ -Rtran V J R> Rtr^„ -^tran V 

The continuity equation (Equation ((5])) becomes 

drN + VDrN = , 

where V = dR/dr is given by equation (|A2|) . 



(Al) 
(A2) 

(A3) 
(A4) 
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Our numerical method is as follows. We choose an initial value for i?tran and an initial form for the function A'^(> R) 
on a logarithmic grid in R, where the grid typically has 400 gridpoin ts p er decade of R. According to Equation (jA4p . 
the value of N remains constant at any R that satisfie s Eq uation (|A2p . Therefore, for each initial N on the grid, 
we integrate its corresponding R according to Equation ()A2|) . We simultaneously evolve i?tran according to Equation 
[3. 

B. SELF-SIMILAR SOLUTIONS 

If one considers either th e co Uisionless or the collisional case, then there is only a single characteristic radius, 
■Rtran- Therefore, Equation (|A4[) may be written in self-similar form when appropriately scaled. We define the scaled 
independent variable to be the ratio of R to its characteristic value, 

X = i?/i?tran(T) , (Bl) 

and the scaled dependent variable to be 

y{x)^Rl,^NiR,T)/r, , (B2) 
with p to be determined below. In scaled variables. Equation (|A4|) becomes 

where 

Self-similarity demands 

^^=^^i?tran-Ar (B5) 

dr 

where ^ is a constant to be determined. Substituting into Equation (jB3[) . we may solve for y to obtain 



y{x) = y(l) <j X > 1 ' ^ ^ 



where we must have A > 1. Equation (|A3p becomes 



1^ = --- -RlZ ( I -'dy + [ x'dy) . (B7) 

^tran S a \Jx<l J x>l / 

It remains to determine the three constants p,A,y{l). To do so, we consider the coUisionless and collisional cases 
separately. 

In the coUisionless case, a/s <^ A/Rtmn- This is equivalent to t he co ndition t ^ ps/{an), i.e., that the time elapsed 
is less than the collision time between small bodies. Equation (|B7p demands that p = 3 since ^ is a constant in 
time, i.e., N/rj = y{x)/ Rf^.^^^. Note that since N is frozen (independent of time) at small x (i.e. for R ^ i?tran), this 
immediately implies that N oc R~^ for R ^ i?tran- Integrating Equation (jB7[) yields a relation between the remaining 
two constants, y{l) and A. We write the resulting self-similar solution as follows: 

^^'^ ^'^^ = -^—9nocoih)f3 i-J^n] , coUisionless, (B8) 

where 

„ , , afl t , 

Rtr.n{t)^———- (B9) 

pa 1 — 3/7 
107^ 

+ 6^ + 277+ 108 



X 



-p 



l-x(l-£) , X<1 
Uix,j)^{ ) p ' ^ . (Bll) 

x-^ (£j , x>l 

In the collisional ca se (a /s » A/i?tran), Equation (jB7|) implies that p — 2, which yields N oc R^^ for R ^ i?tran- 
Integrating Equation (jB7| yields the self-similar solution 

^^^^^ = ^-5coi(7)/2 (1^,7) , collisional (B12) 
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where 



an t 

107(7-3) 

gcoi(7 = 2 I r — TT^ 
7^ + 07 + 16 



and fp{x,'-f) is by Equation (jBlip . 



